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ABSTRACT 

In this article, I discuss W.Kohn's criterion for a metal insulator transition, within 
the framework of a one band Hubbard model. This and related ideas are applied 
to 1-dimensional Hubbard systems, and some intersting miscellaneous results dis- 
cussed. The Jordan Wigner transformation converting the two species of fermions 
to two species of hardcore bosons is performed in detail, and the "extra phases" 
arising from odd-even effects are explicitly derived. Bosons are shown to pre- 
fer zero flux (i.e. diamagnetism) , and the corresponding "happy fluxes" for the 
fermions identified. A curious result following from the interplay between orbital 
diamagnetism and spin polarization is highlighted. A "spin-statistics" like the- 
orem, showing that the anticommutation relations between fermions of opposite 
spin are crucial to obtain the SU(2) invariance is pointed out. 



1 Introduction 



In this article, I will give a brief introduction to the Mott Transition in 
the Hubbard model. This subject is in a state of considerable flux, since it 
is believed by a large number of theorists, to be at the heart of the High 
Tc phenomenon. This point of view was suggested by AndersonS as early 

^ Based on article in Frontiers in Solid State Sciences editors L.C.Gupta and 
M.S.Multani, (World Scientific Publications, Singapore). 
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as in 1986, at the Bangalore ICVF conference. It is less clear as the years 
go by, that the problem has an easy solution. However the validity of the 
Mott-Hubbard point of view as a starting point, in the spirit of providing 
at least the analog of Hq seems very probable to me, of course it may well 
be that there is a mystery Hi remaining to be discovered! The scope of this 
article is much more limited. I would like to describe the transition, and 
discuss some powerful techniques that have been recently brought to bear on 
this problem, within a well defined problem, the 1-d Hubbard model. These 
ideas have also been applied with considerable success in higher dimensions, 
by clever numerical work, and I will direct attention to these articles at the 
end. 

In Section 2 we survey the problem of the Mott transition in various 
Hubbard like models. In Section 3, I discuss the criterion of Kohn for the 
Mott transition, and define the charge stiffness from the Kubo formula point 
of view as well as the energy point of view. In Section 4, I carry out the 
Jordan Wigner transformation and identify extra phases arising from it. The 
diamagnetic inequality for bosons is stated, and an amusing spin statistics like 
theorem for two component hard core boson/fermion system remarked on. 
Amusing odd even effects in diamagnetism and weak spontaneous magnetism 
in the Hubbard model are pointed out. In Section 5, I discuss the Bcthe 
Ansatz equations for the 1-d Hubbard model, which turn out to be very 
useful in extracting the charge stiffness. 

2 The Mott Transition 

We will first of all dispose of the long ranged Coulomb interaction by 
confining ourselves to the single band Hubbard model : there is no good 
reason to do this beyond achieving simplicity. The conducting state, by defi- 
nition, has metallic screening, making the neglect of the long ranged coulomb 
interaction a reasonable approximation. The nature of the insulating state 
obtained within a purely short ranged Hubbard model is surely incorrect with 
regard to the excitations involving promotion to the upper Hubbard band; 
i.e. optical excitations. 

The situation in say the 2-dimensional square lattice is complicated 
by magnetic LRO. At half filling, i.e. one electron per site, the ground state 
is presumably magnetically ordered and insulating. The usual statement of 
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a Mott-Hubbard gap is in terms of an energy gap towards adding an extra 
particle to a half filled band: the existence of a gap of 0(1) distinguishes the 
insulator from a metal. This formulation, applied to the magnetically or- 
dered Antiferromagnet, is confused by the existence of a Spin Density Wave 
gap, and the subtlety of the Mott-Hubbard gap is swamped by the essentially 
trivial SDW gap. We need to consider situations where the band structure 
inhibits the formation of an SDW state: many possibilities arise (a) The tri- 
angular lattice: here the band structure does not have a strong maximum at 
any Q at half filling although a three sublattice ordered AFM state is the 
ultimate large U fate of the model; (b) Random hoppings: this is a doubly 
hard problem with disorder and interactions and no really convincing stud- 
ies exist so far; (c) The Hubbard model on a special lattice in 2-dimensions 
constructed by Shastry and SutherlandBl : (see Figure 1). This is an inter- 
esting two dimensional lattice with a unit cell consisting of four squares of 
the square lattice with diagonal bonds (existing in only two of these squares) 
running NE-SW in one square and NW-SE in the square which shares no 
bonds with the first. We consider a Hubbard model on this lattice with hop- 
pings if: and t\/2a. The band structure gives a semi-metal at half filling, with 
a quadratic touching of the filled valence band, and an empty conduction 
band, for generic a. At large U , the model reduces to a Heisenberg Anti- 
ferromagnetic model with S* = 1/2 on this lattice, with bond strengths say 
J = At'^/U and J' = 2Ja; the really interesting fact is that the Heisenberg 
model can be solved exactly, and has frozen dimer order for the region a > 1. 
For a < 1 the model has not been solved, and approaches the usual nearest 
neighbour model as a — > 0. This is the only case of an exactly solvable 
Heisenberg model with S* = 1/2 in 2-dimensions that I am aware of, and is 
of interest from a Mott-Hubbard point of view, since one expects the semi- 
metal to insulator transition, to be quite representative of the phenomenon, 
uncomplicated by magnetic LRO. 

We confine ourselves in this article to the 1-d Hubbard model, which 
has the great virtue of being solvable exactly by the Bethe Ansatz ™ , and 
also of exhibiting a Mott Transition, uncomplicated by magnetic LRO. The 
model turns out to possess algebraically decaying spin correlations, which are 
by uncertainty principle arguments, the maximal remnants of magnetic LRO 
in 1-dimension. The tool that we bring to bear on this problem, is a powerful 
and beautiful argument due to W. Kohn§. This argument of Kohn, recently 



3 



sharpened and revived by us B after many decades of neglect, turns out 
to be quantitatively implementable in the 1-d Hubbard model. In essence, 
Kohn's argument is that we can distinguish between a metal and an insulator, 
and hence track the Mott transition, by studying the ground state energy 
as a function of the boundary angle. This involves studying the solution of 
the many body problem with twisted boundary conditions, and happily, the 
Bethe Ansatz method admits this complications and remains tractable 0. 



3 Kohn's criterion for the Mott transition 

We first review the arguments of Kohn, specializing to a 1-band d- 
dimensional lattice fermi system, and deduce the implications- which con- 
cern the utility of studying the effect of twisting the boundary conditions 
in a many body system. The second derivative of the ground state energy 
under variation of the flux is related to the zero frequency physics of the 
electrical conductivity. The economy of the Kohn argument is noteworthy: 
we talk of the properties of the conductivity directly, and do not worry about 
one electron gaps. It should be mentioned that several works have appeared 
in literature, trying to use various indirect probes of the Mott transition, 
including the double occupation, binding between holes and doubles, and 
various fictitious fermi-bose fields. The directness of the Kohn argument, in 
contrast, is particularly appealing. 

Consider a d-dimensional hypercubic lattice of linear dimension L, with 
spinless fermions having a nearest neighbour hopping matrix element 't', 
and with arbitrary density dependent interactions that are lattice translation 
invariant and assume periodic bc's. The case of particles with spin is trivially 
obtained from these considerations, in general the different species can be 
twisted in differing degrees, leading to charge as well as a spin stiffness. We 
now introduce a uniform vector potential A^x, which modifies the hopping in 
X directed bonds by the usual Peierls phase factor t — > texp{±i^/L), where 
$ = L{Axe/hc) and the lattice constant ao = 1. Expanding the exponential 
out we find the perturbed Hamiltonian H' = H — ^j^/L — + 

0(<l>^), where jx = 2tYl sinkxC\Ck and = —2t coskxClCk, and H is the 
Hamiltonian for the interacting fermi system. The energy shift of the ground 
state (g.s.) in the presence of the field is Eo($)-^o(0) = D(^^ / L'^-'^ + 0{^^), 
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or equivalently, we define the charge stiffness D: 

The stiffness constant D is given by second order perturbation theory as 

where < T > is the kinetic energy expectation in the g.s. and -E'o(O) = E^. 
We have assumed that < jx > is zero, this is not always true, and some 
examples are given in Sec(4). Higher order (non quadratic) terms in the 
energy shift formula are important when the energy shift is comparable to 
the energy gaps in the spectrum of H. The latter are 0{1/L) in metals and 
so in this case corrections arise when $ is of 0(1/ 1'-'^-^^^). Level crossings 
would occur and perturbation theory would breakdown for $ of 0{7r). The 
region in $ wherein the derivatives are unambiguously defined, therefore 
shrinks with increasing L in two and higher dimensions. 

We next specialize to A^exp{—iujt) leading to an electric field 

Ex = Ax{iuj/c)x, from which the usual linear response formula gives the 
imaginary part of the ac conductivity 



, 1 V- \<0\jxW>\\E,-Eo 



From Eqs(|D and Eq(^ we see that hm^^o"-^ CTxxi^) = ^-D; also 

2 

lim^^ooi^ o'xxi^) = (iT^iid < >. The high frequency behaviour of 
the imaginary part of the conductivity implies for the real part, through the 
usual dispersion relations, the well knownlil /-sumrule: 



3fJe axxioo) doo = — 3-7 < -T > . (4) 

3 an L°- 

More interesting is the small uj behaviour. From standard dispersion theory 
we separate the low frequency part and conclude that 

3?e axx{u^) = —-[D5{huo) + 7^ E I < > - ^o)' - h'uj')] (5) 
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D, the coefficient of 5{huj) , if non-zero, implies free acceleration or infinite d.c. 
conductivity, which is reasonable here since there is no dissipative mechanism 
in the model at T = O''. The coefficient D is essentially the inverse of the 
effective current carrying mass (for free electrons it is vrpe^/m). Therefore 
the /-sumrule is satisfied by the sum of two terms of the same order, the 
stiffness and the "intraband dipole matrix elements". I would like to 
remark that Eq@ was not written down by Kohn, instead he emphasized 
the imaginary part of the conductivity (i.e. Eq(^), and infact discussed 
metals as systems for which the /-sumrule is violated. A ^ight change of 
view point is advocated in the Shastry- Sutherland paper 0, wherein the 
real part of conductivity, Eq(H) is explicitly stated, and it is stressed that 
there is no violation of the /-sumrule, the zero frequency piece, together 
with the finite frequency piece certainly fulfil it. This point of view makes 
the charge stiffness D to appear very much like the stiffness in a superfluid 
or a superconductor. 

This is all very nice, but we still need an algorithm to calculate 
the stiffness! Here Kohn makes the important observation that changing 
boundary conditions is equivalent to applying the $ field, we can absorb the 
Peierls phases by a pseudo gauge- transformation's' S and shift the effect of 
$ into twisted bc's for the wavefunctions : 

\1>(. . .f + Lx . . .) = exp{i^}'^{. . .f . . .). (6) 

Hence, a calculation of the ground state energy with twisted b.c's suffices to 
give D. 

A crucial point (familiar from Landau's fermi liquid theory) is that 
for a Galilean invariant interacting system, an analogous calculation would 
give the coefficient of 6{huj) in Eq(^) unrenormalized by interactions since 
[jx,H] = ( the ffist term in Eq(|5|) becomes the particle density). For 
lattice fermions the operator commutes with the hopping part of H, but 
not with the interaction piece in general and hence for interacting lattice 
fermions there is the possibility that the two terms in Eq(|^) cancel as some 
parameter is varied, signalling a metal insulator transition. The absence of 
Galilean invariance, thus allows the charge carrying effective mass to vary 
with interactions, and in fact to diverge. 

We thus see that Kohn's argument is a means of obtaining a "trans- 
port coefficient" D, from a calculation of the ground state energy. This is 
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less uncomfortable when we recognize that D is not a characteristic of en- 
ergy dissipation, but rather measures the quantity of freely responsive fluid. 
In fact it is obvious that the stiffness so defined is precisely the superfiuid 
stiffness in say a superconductor. Thus Kohn's metal is a superconductor as 
well! True superconductivity may be distinguished from the Mott-Hubbard- 
Kohn metal by introducing disorder. True superconductors are insensitive 
to disorder by the Anderson theorem, whereas the M-H-K metal would be- 
come a dirty metal, with a finite resistivity, and with energy dissipation. 
This brings us to the work of Thouless© , who brought a closely related and 
successful idea to the case of dirty metals, to distinguish between a dirty 
metal and an insulator, i.e. to study the Anderson transition. Thouless's 
ideas are worked out in a purely non- interacting limit, where one examines 
the various one-electron eigenstates of a random potential. The conducting 
states have a finite curvature with boundary angle, but the curvatures of 
successive energy levels are typically opposite in sign (since the dirty metal 
is crudely visualizable as a "periodic" system with a very short period in 
k-space ), whereby a sum over the energy levels would give a thermodynam- 
ically negligible stiffness. The insulating states (i.e. localized states) have 
a vanishing stiffness, typically exponential in system size with a localization 
length setting the scale. A highly nontrivial part of the Thouless argument 
concerns the disorder average of the stiffness for the extended states, and 
relating this to a truly dissipative object, the conductance. This relationship 
has been explicitized in a recent work by Akkermans and Montamboux0. 

It is worth noting that the relationship between D and the ground 
state energy (g.s.e.) is generalizable to finite temperature, with the free 
energy -F($) replacing the g.s.e in Eq(l). This relationship does not seem 
particularly useful for the Mott problem, since at finite temperatures, D is 
expected to vanish exponentially with system size since there exists a ther- 
mal length scale causing an exponential decay of the one particle correlation 
function. 

I would like to mention a beautiful (if somewhat inaccessible ) arti- 
cle on superfluidity by Leggettllll in 1973, which clarifies the role of twisting 
boundary conditions in normal fluids, solids and superfiuids within the con- 
text of rotating bucket like experiments, and illustrates the correspondence 
between these systems and between metals, insulators and superconductors. 
Also the work of Fisher et. al.0 makes precise the connection between the 
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helicity modulus (i.e. (i^F(<l>)/(i$^) and the superfluid stiffness. A recent 
paper by Scalapino, White and Zhang 1^ discusses the various hmits g — *• 
and uj —>■ 0, for the electromagnetic response kernel in 2 and higher dimen- 
sions, with the object of clarifying the role of the shrinking region of $ in 
the Kohn argument alluded to above. 



4 1-d Hubbard model: Generalities and Cu- 
riosities 

In this section we specialize to the 1-d Hubbard model, which was solved 
by Lieb and Wu in 1968. Before discussing the Bethe Ansatz method, 
we discuss some general ideas concerning the stiffness, and the problem of 
particles moving in a flux. We pay particular attention to the Jordan Wigner 
transformation in this problem, which introduces a certain "flux" of its own, 
leading to occasionally surprising results. It is convenient, and natural, to 
convert the model to one involving hard-core bosons. Let us recall the Jordan 
Wigner transformation between Pauli matrices a", with 1 < n < L and 
fermions c^: 

with = 2c\^Cn — 1, and the hermitean conjugate of the Eq(|^. It is straight- 
forward to see that 

ciCn+l = -f^>n+l (8) 

and 

clc = -(-l)^-^+V^ar. (9) 

where N = J^icjci, is the number operator. 

We are interested in converting a fermionic hopping hamiltonian hav- 
ing density dependent interactions, and with a boundary angle (or flux) $j, 
and with cl+i = ci, 

L 

Hf = -J2 (4cn+l exp(?$ f/L) + h.C.) + Vint (10) 
n=l 

to a hardcore boson hopping problem with a boundary angle $6, and cr^+i = 

L 

Hb = - ^(fT+fTn+l exp(i$f,/L) + h.C.) + Vint. (11) 
n=l 
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We can do this in two stages, firstly use the Jordan- Wigner transformation 
Eq(^, and next a unitary transformation Ui = X^n exp(mx/L), with an 
appropriate x, together this yields 

$^ = $j, + 7r(iV- 1); mod27r. (12) 

The reason for preferring to work with bosons is that there are no extra 
phases involved. In fact we can show easily that the boson problem Eq(|TT]), 
generalized to any dimensions, and on any lattice, has a preference to having 
zero flux, i.e. 

> E°(0), (13) 

where E^{^) is the ground state energy in the presence of the flux $. This 
diamagnetic inequality follows from writing the eigenvalue equation for 
Eq(|lT]) in a standard basis, with particles living at various sites, whence the 
phases make their appearance in the off diagonal elements of the matrix. 
We can decrease the ground state energy by increasing the absolute value of 
the off-diagonal matrix elements by virtue of the Rayleigh-Ritz variational 
principle. Inequality (|l^) follows on using | cos(6')| < 1 . Unless the magnetic 
field is a pure gauge field, we must maximize at least one "bond" , and hence 
the Inequality (|13D, is in general, a strict one and may be regarded as a 
statement of diamagnetism of bosons. A corresponding theorem for fermions 
does not exist. In fact the Eq(p!2|) tells us that given the preference of bosons, 
the fermions like (i.e. the minimizing boundary angle is) to have "happy" 
phases: 

= 0, TT for N= odd, even. (14) 

For the Hubbard model, we need a two component Jordan Wigner transfor- 
mation, converting bosons a", r", to fermions c„^|,c„^|, through 

• • • <-i^n = 4,T (15) 
. . . c^^^ . . r„^_ir+ = <^ (16) 

with = 2cl^^-^Cn^i — 1, and = 2c|^ |C„^| — 1. The corresponding phase 
relations are easily seen to be 

= $, + 7r(iV^ - 1). (17) 

$/,i = <l>. + 7r(iVi-l). (18) 
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with obvious notation. We have thus converted the two fermi component 
model to one with two hardcore boson species, with given phases. It may 
appear that the phases of the fermions were completely disposable and not 
very fundamental. However, at this stage, I would like to inject an essentially 
elementary digression into the discussion, showing how important the fermi 
phases really are. 



Digression: Absence of SU(2) for hard core bosons 

The global SU(2) invariance of the Hubbard model is well known, and 
together with the particle hole symmetry, gives the model SO (4) symmetry, 
as recently emphasized by Yang and ZhangM. This represents the (maxi- 
mal) internal symmetry of the system. One may cavalierly expect that the 
hard core bose representation also possesses SU(2) symmetry. This is not 
so. An elementary demonstration of this follows from the examination of the 
commutation relations of the objects: 

= Eicrt - Tn/2 (19) 
1 

= E(a+rr), (20) 
1 

and its hermitean conjugate. The commutator of T" and does not follow 
the usual SU(2) relations, and nor do T" commute with the two component 
boson version of the Hubbard Hamiltonian . The SU(2) relations of course 
follow, if one reinstates the Jordan Wigner strings. This is a kind of "spin- 
statistics" theorem in nonrelativistic physics. 

Getting back to our main theme, we examine the implications of our 
discussions for the stiffness. The fact that the minimizing flux for fermions is 
or TT implies that we need not necessarily expect a positive definite stiffness 
D for fermions, in contrast to the bosonic case. In all cases, we should 
examine the behaviour for an even number of particles closely, since here, 
the state with flux tt is the minimum. The XXZ model, i.e. the Heisenberg 
Ising model0, can be interpreted as an interacting one component fermi gas. 
When the number of particles N is even, one has a non-zero momentum in the 
ground state, i.e < jx >7^ 0. For example think of 2 non interacting particles 
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in a 6-membered ring, we have a degeneracy with momentum ±7r/3, and not 
much can be said about the stiffness in general terms, since the energy shift 
has a hnear term in $. An exphcit calculation using the Bethe Ansatz was 
given for the stiffness in this model in 0, which may be consulted for details. 
The model has an interesting Metal-Insulator transition, with a crystalline 
solid phase, and the stiffness has a jump discontinuity at the transition point, 
very much like the superfluid stiffness jump in the 2-d classical x-y description 
of a bose fluid. 

In the case of the 1-d Hubbard model, the situation is much richer. 
The method of calculating the stiffness will be discussed in the next section, 
here one may just assume that numerical calculations have been performed. 
At half filling, i.e. = Nl = L/2, an interesting phenomenon was first 
noted by Scalapino and co workers 0. They noted that the stiffness D is 
infact negative, for L = 4* integer, but is positive for L = 4* integer + 2. 
The stiffness can be calculated for arbitrary system sizes , and infact, van- 
ishes exponentially with system size with a characteristic behaviour D = 
exp{—L/^), with a known ^, as one would expect for an insulator. An ex- 
planation of the sign was provided using the Aharonov Bohm effect of a 
hole and a double being created, and transported around the Hubbard ring. 
A more direct understanding of this result follows from the considerations 
above. At L = 4 * integer, the number of up and down particles is even, 
and hence the bosons are "unhappy". The state at zero flux is higher in 
energy than that at flux vr, and hence the stiffness is not necessarily positive. 
The fact that it is negative, does not follow from this ( the energy could, 
for instance, have a local minimum at flux 0), and requires actual calcula- 
tion. These considerations, turn out to be not as academic and exotic as the 
reader may suspect. The negative stiffness is infact known, and has been 
observed0 in aromatic organic compounds, which are longer analogs of 
Benzene. The compounds are [16]-annulene and [24]-annulene, where NMR 
shows orbital paramagnetism, and these compounds are essentially the Hub- 
bard model at half filling for L=16,24. 

Another set of curiosities concern the case when the number of par- 
ticles is 4 * integer, away from half filling, in the Hubbard model. At half 
filling, the ground state is a singlet due to a theorem of LiebGHl. Away from 
half filling, however, whenever the number of particles is 4 * integer, the 
ground state is a spin-1 state. Figure 2. shows the singlet and the triplet 
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state energies for a ring of 8 sites, with 8, 6 and 4 particles. The noteworthy 
case of 4 particles has a triplet ground state. For L > 8, similar results hold. 
This result comes about due to an interesting and subtle effect: the spin zero 
state has orbital frustration, i.e. unhappy phases in the bosonic representa- 
tion, which is relieved by partially spin polarizing the system. The ground 
state can be easily shown to have \Stot\ < 1, by an argument analogous to 
that in the Lieb Mattis theoremEyi. We can go to a subspace where the 
Sz = I, having an odd number of up and down electrons, and in this sector, 
there are no frustrating phases, and hence a unique ground state exists. This 
state cannot be orthogonal to the g.s. of the non interacting case since the 
interaction only affects the diagonal matrix elements; for zero U the ground 
state is degenerate between 5* = 0, 1. Hence the spin of the state could be 
zero or unity, but no higher. For infinite U, it is readily shown that 5* = 1, 
from the Bethe equations which become trivial (see Eq(pTD), whereby one 
expects S = 1 for all U. The numerics confirm this. It is fascinating to spec- 
ulate that this spin polarization, arising from orbital paramagnetism, may 
be a more widely useful concept. 



5 1-d Hubbard: Bethe Equations 

We next consider the Bethe Ansatz equations, for the 1-d Hubbard 
model, in the presence of the fluxi, generalizing the Lieb Wu equation o. 
Firstly, it is necessary to show that the Bethe Ansatz goes through in the 
presence of the flux: this is not obvious and should be checked. We checked 
it, and the details are availableB. Basically, the fact that the R matrix in 
the problem conserves the nuber of particles, is sufficient to guarantee the 
Bethe solvability. The equations are, 

M 

Lkn = 2'Kln + $( + 2^arctan[4(Aj - sin(fc„))/f/], (21) 

i=i 

and 

N M 

2 arctan[4(Aj-sin(A;„))/[/] = 27rJ, + (<l>(-<l>()+2^ arctan[2(Aj-Ai)/f/], 

n=l i=l 

(22) 
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where L,N,M are respectively the numbers of sites, particles and downspin 
particles. The quantum numbers 

/„ = integers (1/2 odd integers) : M even(odd); (23) 

Jj = integers(l/2 odd integers) : N — M odd(even). (24) 

The total energy is given by Etot = — 2 cos(/c„). 

The case of $| = — gives us the spin-stiffness in the problem, 
which turns out to be easily computable, and related to the spin susceptibil- 
ity, by a fairly rigorous calculation H. The other case $| = giving the 
charge stiffness, is not computable with the same rigor, since we note that 
the energy change due to the twist is not extensive, nor is it 0(1), it is infact 
0(1/L) in total energy. Morover, subtle level non-crossing effects can and do 
happen when the flux is of 0{7t), and one may not be able to continue 
a solution found at zero flux to a large value of the flux. The usual trick of 
converting the sums to integrals does not guarantee the accuracy of the result 
to 0(1/L), and the answers obtained by this technique, should be verified by 
other means, such as numerical studies. Solutions of the integral equations 
can be found by manipulating the integral equations, assuming that they 
hold for large enough flux. The justification of these manipulations seem 
to require at least integrability (which we happen to have for the Hubbard 
model'0), and also the avoidance of unfavourable fractions in filling. The 
theory of handling these transcendental equations is not qu ite complete at 
this moment, however the careful work of Woynarovich and others has 
made a promising start in this direction. Also solutions using ideas from 
conformal invariance give useful explicit answers . The conclusion is that 
the stiffness vanishes as we approach half filling linearly in the departure 
from half filling. This implies that the Mott transition is directly visible in 
the charge stiffness. 

Similar ideas can be applied in say 2-dimensions, and interesting 
work of Millis and Coppersmithc9 , and Poilblanc and Dagotto© shows 
that the 2-d Hubbard model, as well as a t-J model show similar signs of 
Mott transition. Further applications to a general interacting electron gas 
have appeared recentlyS' . The utility of these calculations in explaining the 
optical conductivity experiments^ in Hi-Tc is noteworthy, and together 
with the the /-sumrule, the charge stiffness gives a good insight into the 
transfer of optical weight observed, within a one band Hubbard description. 
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Figure Captions 



1. The special lattice in two dimensions for which the Heisenberg s — 1/2 
antiferromagnet is exactly solvable. 

2. The spectrum of the Hubbard ring with 8 sites. The numbers (N,S) 
represent the number of particles and the spin of the state. 
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